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Abstract 



We consider localized modes (discrete breathers) of the discrete nonlinear Schrodin- 
ger equation = V'n+i + V'n-i — + cjV'nPV'n) "7 = ±1, n G Z. We study the 
diversity of the steady-state solutions of the form ^n(i) = e*'^*v„ and the intervals 
of the frequency, a;, of their existence. The base for the analysis is provided by 
the anticontinuous limit {oj negative and large enough) where all the solutions can 
be coded by the sequences of three symbols "0" and Using dynamical 

systems approach we show that this coding is valid for uo < u)* k, —3.4533 and 
the point is a point of accumulation of saddle-node bifurcations. Also we study 
other bifurcations of intrinsic localized modes which take place for lv > uj* and give 
the complete table of them for the solutions with codes consisting of less than four 
symbols. 
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1 Introduction 



The discrete nonlinear Schrodinger equation (DNLS) 

^~tt^ = ^2lpn + (7 I r V'n, (l) 

where ipn = 'ipni't), = V'n+i ~2'0„ + '0„_i, and a — ±1, is one of the basic 

lattice models which appeared in last decades in various contexts of physics 
and biology and has been intensively studied (for review see e.g. [1,2,3,4]). 
One of the recent applications of (1) is related to the theory of Bose-Einstein 
condensates in optical lattices, where the mean-field Gross-Pitaevskii equation 
with a periodic potential is reduced to (2) within the framework of the so-called 
tight- binding approximation (see e.g. [5] and references therein) . A particular 
attention has been paid to specific solutions of (1) which are spatially localized 
and periodic in time and called intrinsic localized modes or, alternatively, 
discrete breathers. For numerous applications of model (1) it is desirable to 
have at hands a practical guide where possible forms of intrinsic localized 
modes were described together with the information about their stabihty. The 
aim of this paper is to make a step toward the compilation of this guide. 

The simplest class of the discrete breathers can be obtained using the steady- 
state ansatz ■0n(^) = e"^*iin- Then the stationary amplitudes Un, n = 0, ±1, ±2... 
satisfy the lattice equation 

^2Un + + CT | |^ Mn = (2) 

and the condition of spatial localization 

ii„ — > as n ^ ±oo. (3) 



Eq.(2) is quite general and appears independently in other applications, for 
example, in polaron theory as the equation for stationary localized solutions 
(i.e. polarons) derived on the basis of the semiclassical Holstein model [6]. 
Generically, a can be both positive and negative, uj can be any real value and 
the amplitudes n e Z can be complex. However, for localized solutions (3) 
the analysis can be simplified by taking into account the following points: 

(a) One can show [1] that all the solutions of Eq.(2) which satisfy the local- 
ization condition (3) have the form Un — VnC^^, where n e Z, i>„ and 9 are 
real. 

(b) The localization condition (3) implies that lj does not lie in the phonon 
band defined by < a; < 4. Solutions for a; < and for a; > 4 are connected 
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by the following staggering transform: if discrete function {«„}, n G Z is a 
solution of (2) for a = ctq, cu — cuq < then {(— is a solution of (2) for 
a = —(To and a; = 4 — ujq. 

(c) One can show [7] that for uj < and a — —1 (or, correspondingly, for 
u! > 4 and a = 1) all the solutions of (2) except the trivial zero solution are 
unbounded i.e. | m„ |— > oo as n tends either to +oo or to — oo. 

Thus, for the study of localized solutions Eq.(2) can be replaced by 

AsMn + UJUn + ul = 0, (4) 

where the discrete function n e Z is real and a; < 0. 

Solutions of (4) have been widely discussed in physical and mathematical liter- 
ature. Apart from finding important particular examples, such as the Sievers- 
Takeno mode [8], the Page mode [9], the twisted modes [10,11], some authors 
have addressed the problem of the existence and stability of more complex 
objects [10]. Recently, the analysis of general structure of the set of localized 
solutions of Eq.(4) was proposed in [12,13]. In particular, Ref. [12] is devoted 
to classification of symmetric localized solutions of (4) while in Ref. [13] a 
quite general approach for the study of both symmetric and non-symmetric 
solutions of Eq.(4) was elaborated. Both these papers offer a coding for the 
solutions of Eq.(4). Some discussion of these results in comparison with our 
ones can be found below in Sect. 5. 

In the present paper we study different types of localized modes which can 
coexist at a given parameter u. The general scheme of our approach is as fol- 
lows. For large values of we describe solutions in terms of coding sequences 
using anticontinuous limit approach [14]. To this end we rewrite Eq.(4) in the 
form 

aA2Vn -Vn + v^^O, a= -l/u. (5) 

where = yJ—uiVm n ^ Z. Next, we calculate the boundary cu* until which 
this coding works. Finally, we study bifurcations of the solutions, describing 
them in terms of their coding sequences. This step is implemented using the 
representation of Eq.(4) as a map in M^. In this part our work was influenced 
by the paper [15] where a similar problem was analyzed for the Henon map. 
The information about the bifurcations of the simplest solutions (the solutions 
which have codes of length smaller then 4) is collected in Table 1. 

The paper is organized as follows. In Section 2 we describe the coding of the 
solutions which comes from anticontinuous limit and give an estimation for the 
interval where it is valid. In Section 3 we study the problem from the viewpoint 
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of dynamical systems introducing a 2D map T associated with the discrete 
equation (5). In these terms the problem is reduced to the analysis of the 
intersection points of stable Ws and unstable Wu manifolds of zero hyperbolic 
fixed point of T. We estabhsh the ordering of points of intersection on these 
manifolds. Using this information in Section 4 we describe the homoclinic 
bifurcations. Section 5 contains a summary and discussion of the results. 



2 Anticontinuous limit 

The anticontinuous limit (ACL) corresponds to the zero value of the couphng 
constant a in Eq.(5). When a = the amplitudes of the sites Vn are inde- 
pendent from each other and can acquire one of the three values as follows: 
and ±1. Thus any solution {fn} can be coded by an infinite sequence {sn}, 
n G Z, where Sn is one of the symbols, — , 0, and +, corresponding to Vn = —1, 
Vn — ^ and Vn — i- Localized solutions in this case can be identified with 
those of the sequences which start and end with infinite number of consec- 
utive zeros and contain a finite number of nonzero elements (an example is 
{. . . , 0, 0, +, — , +, 0, 0, . . .}). We denote the set of all these sequences by Cq. 
It is convenient to define the code and the length of the sequences from Cq. 

Definition 2.1 A code of a sequence {rj} e Cq is the maximal subsequence of 
{rj} which has nonzero first and last elements. 

Definition 2.2 A length of a sequence {r}} G Cq is a number of symbols in 
its code. 

Definition 2.3 A word is a sequence of the symbols {—,0, +}. A word which 
consists of a finite number of symbols is called a finite word. 

The above mentioned example {. . . , 0, 0, -|-, — , -|-, 0, 0, . . .} has the code {H — 
-|-} and its length is equal to 3. 

It is known [14] that all solutions which exist in the limit a — can be contin- 
ued by a up to some positive value a*. Moreover, the solutions corresponding 
to the sequences from Cq remain localized under this continuation and their 
decay as n —>■ ±oo is exponential. A basis for the above statements due to [14] 
is given in Appendix A. 

In literature there are various bounds for the limit value a*. Application of 

Theorem 9 of the paper [16] to Eq.(5) yields a* > 1/(10 + 4v^) ^ 0.0639. 
The theorem from Section 3.3 of [1] states that a* > l/6\/5 ~ 0.0745. Below 
we show that a* > (3\/3 — l)/52 pa 0.0807 by making use of the following 
statement 
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Lemma 2.1 For each hounded solution v = {vn\, n of Eq.(5) the fol- 
lowing estimation is valid: 



sup \ Vn — V, 
neZ 



^1 < 4a\/ma. (6) 



Proof. Suppose that v is a solution of (5), (3) and there exists po > 4:a\^l + 4a; 
and no such that |w„(,— t';^^^ | > po- This impHcs that |t'rt(,+i— 2f „(,+fno-i| > Po/ct 
and, consequently, there exists m e {no — 1, no, no + 1} such that \vjn\ > Po/4q! 
and 



Repeating the same procedure one can construct the sequence po, Pi, P2, 
such that 



Pk+l — 




By simple graphical arguments one can show that if po > 4q;^/1 + 4q; then 
Pfc — 00 monotonically as A; ^ 00. Since sup„ |f„ — vf^\ > pk for any k — 
0, 1, 2... we arrive to the contradiction which proves the relation (6) ■. 

Let a be fixed and v = {vn}, n e Z be a bounded solution of (5). Define the 
function 

M{a) = supsupminjlw^l, |f„ — 1|, |f„ + 1|}. (7) 



where the first supremum is taken over all bounded solutions of (5). Lemma 
2.1 implies 

Corollary 2.1 The function M{a) is correctly defined and M{a) — >• when 
a ^ 0. 

Using Lemma 2.1 one can prove also the following statement: 

Theorem 2.1 For some interval < o; < a*, where a* > = (3\/3 — l)/52 
there exists one-to-one correspondence between the set of all localized solutions 
of Eq.(5) and the set of all sequences from Cq. 

The proof can be found in Appendix B. Reformulating this statement for 
initial equation (4) results in the following theorem 
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Fig. 1. Some simplest modes of Eq.(5): (a) the Sievers-Takeno mode (the code {+}); 
(b) the Page mode (the code {++}); (c,d) the twisted modes (the codes {H — } and 

{+0-}). 



Theorem 2.2 For some interval —oo < u < uj* , whereto* >ojq — — 52/(3\/3— 
1) ^ —12.3923 there exists one-to-one correspondence between the set of all 
localized solutions of Eq.(4) and the set of all sequences from Cq. As uj ^ —oo 
the amplitudes of the sites denoted by {±} tend to ±^/—uj, whereas the am- 
plitudes of the sites denoted by {0} tend to 0. 

Theorem 2.1 establishes a correspondence between the localized solutions of 
(5) and the codes of the sequences from jCq for « "close " to anticontinuous 
limit. We call this correspondence ACL coding. As examples of such coding we 
mention (see Fig.l) that the Sievers-Takeno mode has the code {+} (or { — }), 

the Page mode can be coded by {++} (or { }) and the twisted modes of 

two various types have the codes {H — } (or { — h}) and {+0—} (or {—0+}), 
correspondingly. 

The estimation for a* given above is not optimal. In the Section 4 we calculate 
a* numerically. For this we use the dynamical system approach. 
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3 The dynamical system approach 



3.1 Basic notations and some auxiliary statements 



It is convenient to present the real-valued second-order difference equation (5) 
in a form of a two-dimensional area-preserving map T [1,3,17,18] 



T : <^ 



(8) 



X. 



The inverse map T ^ is given by 
: 



x^{2-r'--lx^)x-y-, 



(9) 



y = x. 



In what follows we call a point b = Ta e the T -image of the point a e M^. 
The T-image of a curve in is defined as a set consisting of T-images of all 
points of the curve. 
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The map T is associated with Eq.(5). Specifically, if a = (a;o, yo) ^ is an ar- 
bitrary point of the plane, then the bi-infinite vector formed by x-coordinates 
of points ^"a, n e Z satisfies Eq.(5). The map T admits the following inver- 
sions 



X 



( 



T 



—X 

\-y. 




(10) 



For a > and a < —1/4 the point O = (0, 0) is a hyperbolic fixed point (see 
e.g. [1]) of T. It possesses invariant stable, Ws, and unstable, Wu, manifolds. 
Both Wu and Wg consist of two branches (see Fig. 2): Wg — U W~ and 
Wu = W+ U W~ . We denote by the branch of manifold Wu which goes 
out from the hyperbolic fixed point O keeping y positive. Correspondingly, 
W^ is the branch of manifold Ws which enters the hyperbolic fixed point O 
with positive x. The branches are connected by the symmetries: IV^^ and W~g 
are symmetric with respect to O whereas W^u ^u,s symmetric with 
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Fig. 2. An example of the branches of stable Wg and unstable manifolds of 
the hyperbolic fixed point O. The point H, the intersection of and W^, is a 
homoclinic point. The circles mark the homoclinic orbit associated with it. This 
orbit corresponds to the Page mode for Eq.(5). 

respect to the axis y = x. So, any of the four branches allows one to reconstruct 
three others. Some analytical approximations for W^^ can be constructed 
using normal form analysis, see [18]. 

Points of intersections of Ws and Wu {homoclinic points) correspond to local- 
ized solutions of Eq.(5). As it was shown above, when < a < a* all localized 
solutions can be coded by the sequences from Cq. The homoclinic points also 
can be associated with these sequences. Each homoclinic point has the only 
matching sequence from Cq, but each coding sequence from Co corresponds to 
an infinite number of homoclinic points connected with each other by itera- 
tions of T [homoclinic orbit). This is a natural consequence of the translational 
invariance of the lattice. Due to the involutions (10), (11) the points of inter- 
sections of (and W^) with the lines y = x and y = —x are homoclinic 
points. Their codes are, correspondingly, symmetric or antisymmetric. 

Let us introduce the following definitions. 

Definition 3.1 We call basic points the nine points (0,0), (—1,0), (1,0), 
(0,1), (-1,1), (1,1), (0,-1), (-1,-1), (1,-1) on the plane (x^y). 

Definition 3.2 The spot Ssk, s,k E {—, 0, -|-} is a closed square on the plane 
[x, y) with the side 2M{a), where M{a) is given by Eq.(7), centered at the basic 
point coded with the symbols s and k. The symbol s is the sign of abscissa of 
the basic point (or zero), the second one is the sign of its ordinate (or zero). 

For example, S'+_ is the spot which is centered in the point (1,-1) (see Fig. 3). 
It follows from Corollary 2.1 that 
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Fig. 3. The spots Ssk and stable (dashed) and unstable (solid) manifolds of the fixed 
point O. Small circles mark the basic points. 

(i) for a small enough each homoclinic point on the plane {x, y) is situated in 
the vicinity of one of the basic points. When a tends to zero the distance 
between the homoclinic point and the corresponding basic point also tends 
to zero; 

(ii) all homoclinic points which tend to a basic point as a goes to zero can be 
covered by the spot centered at this basic point. 

(iii) if a is small enough the spots corresponding to different basic points do not 
intersect each other. 

An important fact that follows from (8) is that T maps the homoclinic points 
from a spot Ssk, s,k G {—,0,+} to points belonging to the three spots Sk-, 
Sko, Sk+ but not to any other spots. 

For what follows we need the definitions of segment, initial segment and funda- 
mental segment of W^j^^ (see [15,19]). They can be found in Appendix C. Also 
it is useful to introduce a natural parameterization on the branches and 
W^. Let the parameter r be the length of the curve measured from the 
fixed point O. The point O corresponds to r = 0. The coordinates of the points 
on the branch will be denoted by a;^(r) and y'^ir). In the same manner 
we introduce the natural parameterization on the other branches W~,W^ and 

W-. 

u 

Now we prove the following auxiliary statement: 
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Lemma 3.1 Let 7 = {x{T),y{T)), r e [''■i,T2] be a curve on the plane {x,y) 
such that 

(i) x(r),^(r)eCi[ri,r2]; 

(a) for T e [ti,T2] the curve 7 = (x(T),y(T)) lies entirely within one spot Ssis2) 
si,2 e {-,0,+}; 

(Hi) for T G [ti, T2] the functions x(r) and yir) are monotonic with respect to t; 
(iv) |i/(r)/i(r)| < 1/2 for all r e [ti,T2]. 

Let^ — (x(T),y(r)) be T~^-image of^. Then for a small enough the following 
statements are valid: 

(a) forr G [Ti,r2] the functions x{t) andyir) are monotonic. Moreover, if S2 ~ 
{0} then sign x{t) = sign x{t) and if S2 7^ {0} then sign x{t) = —sign x{t); 

(b) \y{T)/x{T)\ < 1/2 for all r G [ti,T2]. 

Proof. Consider the case when the curve 7 hes entirely in one of the spots 
Ssj^s2 with S2 = {0}; the cases of S2 = {+} and S2 = { — } can be analyzed in 
the same manner. Differentiating (9) one obtains 



x — [2-\ X )x 

a a 



y^x. 



(12) 
(13) 



The side of the spot Sgj^o is equal to 2M{a), so x'^ < 4M^(q;). This implies 
that for a small enough the derivatives x(t) and ?/(r) do not change sign and 
sign x{t) = sign x{t). Dividing Eq.(13) by Eq.(12) and taking into account 
(iv) one obtains (b) ■. 

Consider now the branch W^. Let z = (x^(ro), y^(ro)) be the point of the 
first intersection of with the boundary of the spot 5*00- Then the initial 
segment W^{0, z] remains entirely within the spot 6*00 and on it the following 
Lemma is valid: 

Lemma 3.2 For < t < tq and a small enough the functions x^{r) and 
ytir) CLf^ monotone increasing and |?/^ (T)/a;+(T)| < 1/2. 

Proof. First, we note that in some small neighborhood of the point O the map 
T can be well approximated by its linearization. This means that there exist 
T* < To and the constants Ci^2 > such that for < r < r* 



x+(r) -rcos^l <Cit2 < Cir*^; 
|2/+(r) -rsin99| <C2r2 < C2T*', 
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where is the angle between a stable eigenvector of the linearized map and 
positive direction of the x-axis and 




l + ^)'-l = 0(a), (a^O). 



Let z* = {xf{T*),yf{T*)). Then for a and r* small enough the initial segment 
W^{0,z*] satisfies the conditions of Lemma 3.L The transformation 
maps W^^{0,z*] to another initial segment of for which according to 
the Lemma 3.1 the conditions of the same Lemma 3.1 hold also. Repeating 
the procedure iteratively one concludes that for the whole initial segment 
W^{0, z] the statements of Lemma 3.2 are valid ■. 

Using Lemma 3.1 one can prove the following statements. 

Lemma 3.3 For a small enough there exists a fundamental segment Bf e 
such that 

(i) Bf intersects the spot S+o o-nd does not intersect the other spots; 
(a) the intersection Bf fl Sj^q consists of one connected component; 
(Hi) all homoclinic points which Bf contains belong to S+q; 
(iv) i/h e 5+ is a homoclinic point then for any n T"h e 5*00 • 

Proof. Consider the fundamental segment VF+(T'z, z] where z = (2;(to), ?/(ro)) 
is defined above (see also Fig. 4). Evidently Wg"^(Tz,z] C 5*00 and for it the 
conditions of Lemma 3.1 holds. Define B^ = T^^W^{Tz,z] = l^/(z,T~^z] 
(see Fig. 4). This fundamental segment can be given in a parametric form, 
specifically, B^ = (a;+(r), |/+(r)), r e (To,f] for some f and according to 
Lemma 3.1 for tq < r < f the coordinate xf{T) grows monotonically and 
\y'^ (t) / {t)\ < 1/2. Then using simple geometrical arguments one arrives to 
(i) and (ii). The point (iii) is a consequence of the fact that all homoclinic 
points belong to one of the spots. The point (iii) implies also (iv) ■. 

By symmetry one can conclude that for a small enough there exists a funda- 
mental segment B~ C W~ such that all homoclinic points which belong to 
B~ lie in S^q. Analogously, there exist the fundamental segments B^ C 
and B~ C W~, such that all homoclinic points which belong to B^ lie in So±. 
We call Bf and B^ the basic fundamental segments. 

The following statements are valid 

Lemma 3.4 For a small enough and for any n e Z there exist the following 
one-to-one correspondences: 

(i) between the homoclinic points which belong to the fundamental segment 
T'^Bf and the set of the sequences from Co which have the codes with last 
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{-+} 



'{0+} 











Sof 





















r {++} 

2 



T" 




Fig. 4. A basic fundamental segment = W^{z,T ^z]. All homoclinic points 
belonging to B'^ also belong to r\^^ = 5+o H Bf. T^^-images of these points lie in 
r|°+> and r^++> 



r 



. The notations are defined in the text. 



symbol "+"; 

(a) between the homoclinic points which belong to the fundamental segment 
T'^B~ and the set of the sequences from Cq which have the codes with last 
symbol " 

(Hi) between the homoclinic points which belong to the fundamental segment 
T^B^ and the set of the sequences from Co which have the codes with first 
symbol "+"; 

(iv) between the homoclinic points which belong to the fundamental segment 
T^B^ and the set of the sequences from Cq which have the codes with first 
symbol "— 

Proof. As an illustration let us consider the point (i) and n = (the basic 
fundamental segment B'^ itself); all other cases can be treated similarly. Let h 
be a homoclinic point from B^. Consider the corresponding solution of Eq.(5). 
Its coding sequence from Cq contains the infinite block (+,0,0,0...), so the 
code end up with the symbol "+". From the other side, consider a sequence 
L E Cq which end up with the block (+, 0, 0, 0...). If a is small enough, there 
exists a solution of Eq.(5) which corresponds to L. On the plane {x,y) its 
image is a sequence of homoclinic points hk, k G Z, Thu = hk+i- According to 
the coding sequence, there exists a number fco, such that for k > ko hk & Sqq 
and hkg G 5*0+. Then /i^,, G B^ by construction of B^ in Lemma 3.3 ■. 

In what follows we will label homoclinic points by their codes saying, for 
example, "homoclinic point { — I — }". This labeling is well-defined within 
one fundamental segment: two points which belong to the same fundamental 
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segment definitely have different codes. In what follows if the fundamental 
segment is not specified, this means that the homoclinic point under consid- 
eration belongs to Bf. 



3.2 Ordering on B 



We say that a >s b (a <s b) on the branch if b lies closer (farther) on 
to the hyperbolic fixed point O than a. The ordering on is defined 
in the same manner. 

We denote by {{Ci}{C2}} the word composed by two words {^1} and {^2}- For 
example, if {6} = {++} and {6} = {— } then {{6}{6}} = {+ + — }■ 

Then the following Theorem establishes the ordering on B^: 

Theorem 3.1 Let a be small enough and be a word of arbitrary but finite 
length. Then if the total number of symbols {+} and { — } in {^} is odd then 
for any three words {^i},{^2};{%} the homoclinic points on B^ are ordered 
as follows 

{{vi} - {0+} <s {{mmy+y <. {{^3} + {o+y (14) 



and if {C} contains even number of the symbols {+} and {—} then 

{{Vi} + {0+} <. {{V2}0{^}+} <s {{%} - (15) 

Proof. First, let us verify this statement for the words {^} of length 1. Define 
Tq^^ = Wj'{0, z] (where z is defined in subsection 3.1) and rj;''''^ = B^ fl S^o 
(see Fig. 4). Evidently, TFj;^-'^ C Fg^^ According to Lemma 3.3 r\~^^ consists 
of a single connected component and there exist and such that t\^^ — 



(X 



s 



T),y^{T)) for Ti < r < T2 . Consider now T ^-image of F} ' . According 
to Lemma 3.1 T-ir{+^ = (a;+(r), y+(r)) where r varies within some interval 
fl < T < f2, the function x+(t) monotonically decreases within this interval 
and 

\ytir)/xt{T)\< 1/2, Tt<T<Tt. (16) 



The homoclinic points which lie in rj;^^ are T-images of homoclinic points 
which belong to fundamental segment T~^Bf. They are situated in the spots 
S++, 5*0+ and So, T'^Fi^^ intersects these three spots. Since xf decreases 
monotonically on T''^r{'^^ the intersections are in the order given above (see 
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Fig. 4). It follows from (16) that T-^ri^^ does not intersect any other spot. 
Denote 



rj++>=r-iri+>n,s++, 

rr>=r-M^>n5o+, 
rj-+>=r-iri+^n,s_+. 

Each of r2^^\ ^2^^^ ^^'^ ^2 consists of one connected component. On 
the fundamental segment T^^Bf the arc r2^^-'^ includes all homoclinic points 
corresponding to the codes {{^i} + +}; the arc Fg^^^ includes all homoclinic 
points corresponding to the codes {{^2)0+} and the arc r^'^^ - all homoclinic 
points corresponding to the codes {{773} — h} (here rji, 772 and 773 are arbitrary 
words). This implies that on 



{{Vl} + +} <s {{V2}0+} <s {{Vs} - +}■ 
In general case we denote the corresponding segments by recurrence relations 

r{-Hc}+} = r-^F?jP+>n5_fe, 

where {(} runs all the words of length (n — 2) and k e {—,0,+}. Arcs 

r{+fc{C}+}, rl^HC}+} and F^-'={^>+> belong to fundamental segment T-^^-'^^B+. 
Similarly to the considered case using Lemma 3.1 one can conclude that 

(a) T'^fIj'l^P^^ intersects the three spots S+k, Sok and S^k and does not 
intersect the other spots; 

(b) Each of F^+*^^^^+^, F^°'^^^^+^ and Fj^~''^'^^+^ consists of only one connected 
component, so each of these arcs can be parameterized in the form (x+(t), ^^(t)) 
at appropriate interval of r. 

(c) The order in which F^+'^{^>+>, Tjl^''iO+} and F{-'=^^>+> appear on T-("-i)5+ 
coincides with the order in which F^^^^'^'*'^, F^^P''"-'^ and F^Z^''-'^'''-'^ appear on 
rp-{n-2)-Q+ if ^ ^ |Q| -g opposite to it if A; = {+} or = {-}. 

A schematic picture of relations between the arcs Fj^^^J'+J' is shown in the 
Fig. 5. Taking into account that the homoclinic points which belong to Fj^^^J'+J' 
correspond to the codes {{77}{C}+} where rj is an arbitrary word we obtain 
the statement of the Theorem ■. 
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To 

Fig. 5. A schematic picture of relations between the arcs F, 



Remark. Theorem 3.1 gives a simple recurrence algorithm allowing one to 
establish the proper ordering of homoclinic points with codes of lengths smaller 
than some prescribed length on B^. It consists in the following steps: 

Step 1. To write the starting codes {++}-{0+},{ — h} (exactly in this order) 
and the arrow pointed from { — h} to {++}■ 

Step n. Let the codes for (n — l)-st step be just written. 

Then one has 

1. To draw arrows corresponding to each of the codes such that: 

(a) the extreme left arrow is pointed in opposite direction to the extreme left 
arrow on the previous step; 

(b) the neighboring arrows have opposite directions. 

2. To generate new string of codes replacing {^} by the three codes {— {C}}' 
{0{^}}, {+{C}} if the arrow corresponding to {^} is pointed to the right and 
{+{C}},{0{C}},{-{^}} if the arrow is pointed to the left (see Fig.6). 

In Fig. 6 the largest arrow corresponds to the first step of the algorithm, three 
smaller ones correspond to the Step 2, and nine short arrows correspond to the 
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• ® • • ® • • • — 

{-++) {++} {+++} {+0+) {0+} {-0+} {— +) {-+} {+-+} 



Fig. 6. How to place the homoclinic points on Bf in the proper order (illustration 
to the algorithm). 

third step. As a result, the picture shows the mutual position of the homoclinic 
points with the code length smaller than 3 on Bf. 

The same algorithm allows one to order homoclinic points on other basic fun- 
damental segments Bj , B^ and B~ . By comparing orders of homoclinic points 
on both stable and unstable manifolds, it is possible to represent schemati- 
cally how and W~ intersect W^^ (see Fig. 7). Also it is possible to describe 
a structure of homoclinic points which are situated in a neighborhood of a 
homoclinic point {^} (see Fig. 8). 



4 Homoclinic bifurcations 

Theorem 2.1 states that for the interval < o; < a* all localized solutions 
of Eq. (5) are in one-to-one correspondence with the set of all sequences from 
jCq. li a exceeds the critical value a* this does not take place any more and 
solutions can appear or disappear as a result of some bifurcations. Since the 
localized solutions correspond to the points of intersections of the stable Wg 
and the unstable manifolds of the fixed point O, these bifurcations can 
be found numerically by observing how the forms of Wg and change as a 
grows. In such a way, each of the solutions of Eq. (5) which exists for < a < a* 
can be continued until some value of a > a* where it dies in some bifurcation. 
It is convenient to label this solution with the same code as it has in the 
anticontinuous limit. In this sense, the coding of solutions used for < a < a* 
can be used for a > a* also. A problem can occur, however, if some new 
solutions appear since they have no counterparts in the set of the coding 
sequences. 

In order to study possible bifurcations we used a numerical procedure which 
allows one to construct initial segments of the manifolds Wu and Wg, to find 
their points of intersection, and to identify them with possible codes from jCq. 
Using this program we observed that the deformation of the manifolds Wu 
and Wg as a grows is quite monotonic. This results in numerous bifurcations 
of death of homoclinic points. However, we did not observe any bifurcation 
of birth of new homoclinic point. This allows us to formulate the following 
conjecture. 
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Fig. 7. Schematic representation of the intersections of (thin hne) and W, 
(bold hne) with the branch W^. 
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Fig. 8. The homoclinic points in a neighborhood of the homocUnic point with the 
code {^}. 

Conjecture 4.1 As a grows no bifurcations of birth of new solutions occur. 

This conjecture is similar to "No- bubbles-conjecture" of [15] for Henon map. 
In view of this conjecture the problem is reduced to the description of the 
bifurcations of death of solutions which can be identified with their codes in 
the anticontinuous limit. 

Let us introduce the following notations. Let {^} be a finite word. We denote 
by {^} the word formed by the symbols of {^} taken in inverse order. Also we 
denote {^} the word obtained from {^} by replacing the symbols by " — " 
and vice versa. We say that {^} is symmetric if {C} = {C} antisymmetric 
if {^} = {C,}. Also we call symmetric (antisymmetric) the homoclinic point 
coded with symmetric (antisymmetric) word {^}. In the same manner we 
define symmetric (antisymmetric) solutions of Eq.(5). Then the two types of 
homoclinic bifurcations are generic in our case: 

Saddle-node bifurcation. In a point of the saddle- node bifurcation a = asN 
two homoclinic points which exist for a < asN, merge at a = asN, and do 
not exist for a > asN- It is clear that if two homoclinic points with codes {.^i} 
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Fig. 9. The stable (dashed line) and the unstable (solid line) manifolds of the fixed 
point O. (a) a ~ 0.2632, (u = 3.8). In this case a < a* and the structure of the in- 
tersections of the stable and unstable manifolds is the same as in the anticontinuous 
limit, (b) a 0.3125, (w = 3.2), some homoclinic points disappeared; (c) a 0.5, 
{oj = 2); (d) a w 1.25, {lj = 0.8). 



and {^2} undergo the saddle- node bifurcation, then the pairs of homoclinic 

points ({6)5 {6}), ({^i,^}), and ({^1}, {6}) also undergo the saddle-node 
bifurcation at the same a = asN- 

Pitchfork bifurcation. In a point of the pitchfork bifurcation a = ap two 
homoclinic points with codes {^1} and {^2} related by the symmetry {^1} = 
{^2} or antisymmetry {^1} = {^2} die by merging with a third homoclinic 
point {^3}. The code {.^3} of the resulting homoclinic point is symmetric if 
{^1} = {^2} and antisymmetric if {^1} = {^2}- This bifurcation is generic 
since the dynamical system possesses involutions. 
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Fig. 10. Saddle-node (a) and pitchfork (b) bifurcations illustrated in terms of inter- 
sections of lobes of stable and unstable manifolds. 

The both types of bifurcations can be illustrated schematically considering a 
generic intersection of the lobes of Wg and Wu (see Fig. 10). It follows from 
Fig. 10 that a pitchfork bifurcation is typically followed by a saddle- node bifur- 
cation. In this saddle-node bifurcation the symmetric (antisymmetric) homo- 
clinic point (the "output" of the pichfork bifurcation) merges with one more 
homoclinic point of the same symmetry. We checked numerically (see below) 
that for a > a* the both types of the bifurcations take place. 

Following [15] we call homoclinic points a and b double neighbors if there 
are no other homoclinic points on the segments iy„(a, b) and Ws{sL,h). Ho- 
moclinic points cannot undergo bifurcation unless they are double neighbors. 
This observation implies that 

Lemma 4.1 The codes of bifurcating homoclinic points (solutions) must start 
and end with the same symbol. 

Proof. This statement immediately follows from the fact that the homoclinic 
points with codes starting with " + " lie on the branch whereas ones 
with codes starting with " — " belong to W~ . So, the homoclinic points with 
different first symbol of the code belong to different branches and cannot be 
neighbors. The same arguments can be used for the last symbol of the code 
■. 

Let us now turn to the analysis, what happens on the boundary of the interval 
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< q; < a*. 



Lemma 4.2 The value a = a* is not a value of saddle-node or pitchfork 
bifurcation. 

Proof. It follows from Theorem 3.1 that for < a < a* there exist infinitely 
many homoclinic points between any pair of homoclinic points which belongs 
to Bf . This implies that there are no double neighbors at Bf for a — a* M. 

Let {^} be arbitrary finite word. Then we define the word 



{er ^ mm ■ ■ ■ m 

^ V ' 

n times 
and its limit case 



{er^{---{e}{e}---{e}---}- 

The results of numerical investigation allow us to formulate the following state- 
ment: 

Statement 4.1 The value a* ~ 0.28958 is a limit value of a sequence of 
points of homoclinic saddle-node bifurcations. Bifurcating pairs of homoclinic 
points have the codes {+{ — h}" + {H — }"} and {+{ — |-}"'0{H — }"■}, n - pos- 
itive integer. This value a* also can be understood (and found numerically 
in practice) as a point of saddle-node bifurcation of the solution of (5) coded 
{{ — 1"}°° + {"I — }°°} which is double asymptotic to the periodic solution 
{+-}-. 

We comment Statement 4.1 graphically as follows. In Fig. 11 there are plotted 
the parts of stable and unstable manifolds which pass through the area of 
the bifurcations. There are two rectangles marked as (A) and (B) in Fig. 11. 
The coding of homoclinic points within these rectangles can be understood 
from Fig.8. Specifically, in the case of the rectangle (A) one should take ^ = 
{++} and reflect Fig. 8 with respect to central vertical line; in the case of 
the rectangle (B) one should take ^ = {+0} and rotate the Fig. 8 by 180 
degrees. As a grows the points from the right lower corner of (A) tend to the 
corresponding points from right upper corner of (B). Up to the lines shown 
in Fig. 8 this corresponds to the saddle-node bifurcation of homoclinic points 

with codes {-I h + H — } and {-I hO H — }. Completing the picture of 

Fig. 8 i.e. adding more horizontal and vertical fines we found other codes of 
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Fig. 11. The parts of the stable and unstable manifolds which pass through the area 
of the bifurcations, a ~ a*. The bolder lines correspond to the earlier passages of 
the stable and unstable manifolds through the area. 

bifurcating homoclinic points: 

{- + - + + + -+} 

{+- + - + + + - +-} 
{- + - + - + + + - + -+} 



{- + - + + -+}, 
{+- + -+0 + - + -}, 
{- + - + - + + - + -+}, 



(17) 



Bifurcations of words with greater lengths occur at smaller values of a than 
the ones of words with smaller lengths of code. The codes of the bifurcating 
homoclinic points approach the infinite sequences 

+ and {{-+}°°0{+-}°°}. (18) 



These sequences correspond to the orbits homoclinic to periodic orbit with the 
code {H — }°°. The bifurcations of these orbits can be found by studying stable 
and unstable manifolds of the fixed point {+}°° of the squared map T^. This 
analysis shows that (18) undergo saddle-node bifurcation at a* ~ 0.28958. 
Also it was confirmed numerically that the points of the bifurcations of pairs 
(17) approach this value as the length of the corresponding codes grows. 
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Fig. 12. The saddle-node bifurcation of the solutions with the codes {— +} and 
{ !-+}• The point of the bifurcation is asN ~ 0.772 (cj w —1.295). 

Let us now consider in more detail the bifurcations of these solutions which 
occur for a > a*. The complete description of the bifurcations for the solutions 
with codes of length smaller than 4 are shown in Table 1 which was compiled 
studying directly the transformation of the stable and unstable manifolds Wg 
and Wu as a grows. In this Table we show also the bifurcations related to 
them, for example, pitchfork bifurcations which involves also solutions with 
codes of length 5 and 6. It worth noting that the bifurcating solutions may 
have different lengths of the codes. 

In order to study the bifurcations of solutions with greater lengths of the code 
we used the following approach. First, using the Newton method it is possible 
to follow numerically the continuation by a of any localized solution of Eq. (5) 
coded with {^} from the anticontinuous limit. This can be done until some 
final point a = a^^y where the determinant of Newton matrix vanishes. Let 
at the point a = the amplitudes of the sites obtained by continuation of 
the solution with code {^} be vj^^^{a{^y), n e Z. If these final amplitudes are 
close to each other for two different codes {^i} and {^2} and at the same time 
q;{^j is close to ^{ga} then one can expect that the solutions with these codes 
bifurcate. 

We fulfilled this continuation procedure for the solutions of Eq.(5) which have 
the codes of lengths equal or smaller than 9. Owing to the symmetries we 
restrict ourselves by the codes which end up by " + ", i.e. 3^ = 6561 codes. 
Using these continuation results we created a database. This database includes 
the codes {^} of the solution, the values of and the values for 
large enough number of sites (basically we took 40 sites). The search within 
this database allows by an introduced code of solution to find a code of another 
solution which bifurcates with it. The criterion for codes {^1} and {^2} to be 
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a type codes 

{ +}4»{-0-+}; {+ + +-}4»{+0 + -}; 

0.324 SN 

{- + ++}4»{- + 0+}; {+ }4»{+-0-}; 

0.374 SN {+ + +}^{+o+}; { }^{-0-} 

{+0 + +} ^ {+ + ++} ^ {+ + 0+}; 
0.485 P L J L J L J> 

{_0--}^{ }^{__0-}; 

0.609 SN {+ + ++} ^ {+00+}; { } ^ {-00-}; 



0.745 P 



{+ + - + -}=^{+-+-}^ {+- + --}; 

{-- + -+} ^ {- + -+} ^ {- + - + +}; 

0.745 + A SN {+ + _ + __}^ {+_+_};{__ + _++} ^ {_ + _+}; 



0.747 P 



{+ + - + } ^ {+ _ +} ^ {+ _ ++}; 
{--+-}^{- + -}^{- + — }; 
0.747 + A SN {+_+}^{+ + _ + +}; {_ + _}^{__ + __}; 



0.762 SN 



(+0 -+}^ 1+0 - ++}; {-0 + -}^ 1-0 + }; 

{+ - 0+} + -0+}; {- + 0-} +0-} 



{+ + -}^{+-}^ {+--}; 

0.772 P 

{--+}^{-+}^ {- + +}; 

0.772 + A SN {_+} ^ {__++}; {+_}^ {+ + __} 



0.779 P 



{+ + --+} ^ {+ + --++} ^ {+ - - + +}; 

{-- + + -}^{-- + + —}^ {- + + --}; 
0.779 + A SN {+__+} 4^ {+ + __++};{_ + +_} 4^ {__ + + __}; 



1.091 P 



{+ + 0-}^{+0-}^{+0--}; 
{--0+}^{-0+}^{-0 + +}; 
1.091 + A SN {+0-}^{+ + 0--}; {-0+}^{--0 + +}; 

{+ + 00-} =^ {+00-} ^ {+00 - -}; 

1.290 P 

{- - 00+} =^ {-00+} ^ {-00 + +}; 
1.290 + A SN {+00-}^{+ + 00--}; {-00+}^{--00 + +}. 



Table 1 

Bifurcations of solutions with the code length smaller than 4 and related to them. 

By the symbol A we denote the difference between the values of a for pitchfork and 
subsequent saddle-node bifurcations when they are too close to be distinguished by 
the numerics. 
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Fig. 13. A pitchfork bifurcation. Two asymmetric solution with the codes {+0 + +} 
and {+ + 0+} merge with the solution with the code {+ + ++}. The point of the 
bifurcation is ap ^ 0.485 [cu w -2.062). 

the codes of bifurcating solutions was taken as follows 



Here e is some small positive number (we took e = 0.1). 

We test this method using the values of bifurcation points given in Table 1 
and in all the cases the discrepancy was about 10~'^ which corresponds to the 
step in the continuation procedure. In general, the following situations can be 
distinguished: 

- If a solution dies in a saddle-node bifurcation which is not the second bifur- 
cation in the pair "pitchfork - saddle-node" the search in the database gives 
one code for the solution which bifurcate with it. 

- If a solution dies in a pitchfork bifurcation the search in the database typi- 
cally gives 3 corresponding codes. The reason is that if the initial segment 
of stable (unstable) manifold is large enough, its lobes become very narrow, 
so the points of pitchfork bifurcation can be hardly distinguished from the 
points of subsequent saddle-node bifurcation (in the Table 1 we denote the 
distance between these close points of bifurcations by the symbol A; we 
should stress that A is not a fixed value but any value which is below the 
accuracy of the method). So, the search in the database in this case gives 
two codes of solutions which participate in the pitchfork bifurcation and one 
more code for the solution which dies in subsequent saddle-node bifurcation. 

- By the same reason, if a solution dies in a saddle-node bifurcation which is 
the second bifurcation in the pair "pitchfork - saddle-node" the search in 
the database typically gives 3 codes. 



mm 

|m|<10 



(E (4-m("te}) -^n^'H"te}))' < e. 



(19) 
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Fig. 14. Number of coexisting solutions with codes smaller or equal to 9, N{a), 
versus a. 

Also the database allows to calculate the number N{a) of solutions with codes 
of length smaller or equal to 9 which coexist for a given value of a. The plot 
of N{a) versus a is shown in Fig. 14. The following comments about the 
bifurcations for a > a* are in order. 

(a) There exist intervals of a where no bifurcations occur. One of them is 
marked in Fig. 14 with ellipse. Its boundaries are ai ~ 0.295 and 02 ~ 0.314. 

(b) On the interval {a*,ai) the first sequence of bifurcations occurs. At this 
interval there die all solutions whose codes contain the sequences of symbols 

-+0 + -, - + + + -, +-0-+, + +. (20) 

This fact is confirmed by the analysis of intersections of stable and unstable 
manifolds Ws and Wu- The typical bifurcations here are 

{{^1} - + + + - {V2}} ^ {{vi} - +0 + -te}}; 

{{^1} + + {^2}} ^ {{r/i} + -0 - +{V2}}. 

The end value ai corresponds to the bifurcations 

- + + + - ^ - +o + -{-+}~}; 
{{-+}°° + — + {+-}°°} ^ {{-+}°° + -0 - 
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of the orbits homoclinic to periodic orbit {H — }°°. This situation is similar to 
one found for Henon map [20]. Following [20] we suppose that for ai < a < a2 
the localized solutions of Eq. (5) can have any code except ones which includes 
the words (20). 

(c) One more sequence of bifurcations starts with the value 0:2 which corre- 
sponds to the bifurcations 



{{-+ro+o{-+}°°} 

{{+-}^0 + 0{+-}°°} 

{{+-ro-o{+-}°°} 

{{-+}°°0-0{-+}°"} 



{{-+r + +0{-+}~}; 
{{+-rO + +{+-}°"}; 

{{+-r--o{+-r}; 

{{-+}°°0--{-+}~}. 



These bifurcations affect the solutions which codes contain the sequence of 
symbols " — hO + O" and other related to it due to symmetries. At the interval 
(a*, 0.341) all these solutions die. The typical bifurcation is 



{{vi} - + + ^ {{m} - +0 + o{772}} 



and other related to it by symmetries. 

(d) For a > 1.21 all the solutions in our database have the codes which are 

composed by the words " + " , " — "," + +" , " " separated by groups of the 

symbols "0" . So, they can be regarded as bound states of Page and Sievers- 
Takeno modes. A study of such multipulse solutions for a large can be found 
in [21]. 

Since u = —1/a all these results can be rewritten in terms of Eq.(4). In 
particular one obtains that ACL coding holds for a; < a;* ~ —3.4533 and the 
point cu* is the point of accumulation of saddle-node bifurcations of solutions 
of Eq.(4). 



5 Conclusion 



To summarize, the main results of the paper are: 

(a) For the equation Eq.(4) we used the coding of all localized solutions (3) 
by the words which consist of the symbols "— ", "0" and "+". This coding 
comes from anticontinuous limit (ACL-coding). We calculate the boundary of 
the interval of parameter 3.4533 for which this coding is vahd. 



27 



(b) In order to do this we analysed the problem from the viewpoint of dynami- 
cal systems. Since the localized solutions of this discrete equation corresponds 
to homoclinic points of zero fixed hyperbolic point O of associated map T we 
study the mutual position of stable and unstable manifolds of O. A simple 
rule was found for the ordering of the homoclinic points on these manifolds 
when anticontinuous approximation is valid. 

(c) Based on this we studied the bifurcations of localized solutions of Eq.(4). 
They correspond to homoclinic bifurcations for the map T. It was shown that 
the end point u* of the interval where the ACL-coding is applicable is the point 
of accumulation of saddle-node homoclinic bifurcations. Then we studied the 
bifurcations which occur for uj > uj*. We found an interval of uj where no 
bifurcations occur and gave a complete list of bifurcations for the solutions 
whose coding sequences have lengths smaller or equal to 4. 

It worth to compare our results with the results of [12] and [13]. In both these 
papers Eq.(4) is rewritten in the form 

A„+i + A„_i + C^„ = 3A^. (21) 

The authors of [12] calculated numerically the values of C for which there 
exist symmetric localized solutions keeping fixed two internal parameter cor- 
responding to the asymptotic of the solution for n large (specifically, e and 
N) . The values of these parameters were chosen for reasons of numerical con- 
venience. The coding of solutions introduced in [12], in fact, corresponds to 
ACL coding. Applying our results one can found that for C > C* ^ 5.4530 all 
solutions of (21) can be coded by Cq sequences. However, a// values of C given 
in resulting Table 1 of [12] exceed this threshold. This explains why the num- 
ber of found solutions is equal to one predicted by the analysis of number of 
possible coding sequences (which, generically, does not take place). Changing 
the internal parameters of the method one can make the value of C smaller 
than the critical value C* and the situation immediately becomes much more 
complex. We note also that the ilhistrativc example of [13] related to Eq.(21) 
corresponds to C = 8 which also exceeds the threshold value C*. 

To complete the picture, the information about the existence of discrete brea- 
thers should be accompanied with the results on their stability. This part 
of the work will be presented elsewhere. Here we note that the analysis of 
bifurcations of discrete breathers proved to be quite useful for the study of 
their stabihty. We sketch the relation between these topics as follows. In the 
linear approximation the stability problem can be reduced to the eigenvalue 
problem for non-symmetric infinite matrix L. This matrix is the product of two 
infinite symmetric matrices, L_ and L+, where L+ correspond to linearization 
of Eq.(4) (Newton matrix) [10]. An information about spectral properties of 
L_ and L_|_ can be obtained in the anticontinuous limit. Varying cu one observes 
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that at the point of homochnic bifurcation u = Cj (and not in other points) 
one of the eigenvahies of Lj^ vanishes. From the other side one can show that 
the number of positive eigenvalues of L_ does not change when uj varies. These 
facts, together with the relation between the numbers of positive eigenvalues 
of L_, L+ and L [22,23] give a basis for the analysis of stability of discrete 
breathers. 
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A The continuation from anticontinuous limit 



Let us denote 



(A.l) 



Then 



da 



(9F(v,a) 



(9F(v,a) 



da 



(A.2) 



where 



aF(v,a) 
da 



= col (. . . , A2V„_i, A2V„, AaVn+i, . . .) 



and 



av 



a {-2a-l^?,vl_^) 



a 



a 



(-2q;- 1 + 3^;^) a 



\ 



(A.3) 



are an infinite column- vector and a three-diagonal infinite matrix correspond- 
ingly. The continuation by the parameter a is possible while the operator DyF 
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remains invertible. According to [14] this takes place until some value a — a* . 
The paper [14] offers also a method how to estimate the boundary a* . 



B Proof of Theorem 2.1 

First, we prove the following lemma: 

Lemma B.l For < a. < a^, where ao = (3\/3— 1)/52 fa 0.0807 the operator 
DyF (see Appendix A) is invertible. 

Proof of Lemma B.l. Let us select the conditions when the matrix (A. 3) has 
main diagonal domination (MDD). By definition, MDD means that the mod- 
ulus of each element on the main diagonal is greater than the sum of module 
of other elements of the same string i.e. 

Sn=\-2a-l + 3vl\ -2a>0 

for any n G Z. Having MDD situation one can conclude by the generalized 
Gershgorin theorem that the matrix (A. 3) is non-degenerated. Consider the 
amplitude of one of the sites, v, and suppose that v >0 (the case v <0 can be 
analyzed by symmetry). The value of v can fall into one of the three intervals: 

(i) < t; < l/\/3; 

(ii) l/VS <v<l; 

(iii) v>l. 

The case (i). According to Lemma 2.1 

v-v^ < 4a VI + 4a. (B.l) 

The MDD condition is 



One can check that the condition (B.2) holds automatically if < a < 1/12. 

The case (ii). Let us denote e = 1 — v; then < e < 1 — l/y/S. The condition 
(6) takes the form 



5 = l-3v'^>0. 



(B.2) 



- 3e^ + 2e < AaVTT^. 



(B.3) 
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Prom the other side one has 



S=\-2a-l + 3t;^| - 2q; > St;^ - 1 - 4q; = 3e^ - 6e + 2 - 4q; 
and the MDD takes place if 

3e^ - 6e + 2 - 4q; > 0. (B.4) 

One can check that Eq.(B.4) holds automatically if Eq.(B.3) is satisfied and 
< q; < tto; here ckq = (3-\/3 — l)/52 can be found as a solution of the system 
of equations 

el — Sel + 2eo = 4q;oV1 + 4q;o; 
Sel - 6eo + 2 - Aao = 0. 

The simplest way to make sure of the last statement is to plot in the plane 
{e,a) the areas corresponding to Eq.(B.3) and Eq.(B.4). 

The case (iii). In this case we denote e = 1 + v; then e > 0. The condition (6) 
takes the form 

+ 3e^ + 2e < 4q;V1 + 4q; (B.5) 

and MDD condition is 

5 > 3^2 - 1 - 4Qi = 3e2 + 6e + 2 - 4a > 0. (B.6) 

One can check that Eq.(B.5) imphes Eq.(B.6) for any a > and e > 0. Taking 
the smallest interval between those found in (i)-(iii) we arrive at the result of 
Theorem B.l ■. 

The proof of Theorem 2.1 follows immediately from Lemma B.l and Theorem 
3.3 of [1]. ■ 



C Some definitions from dynamical systems theory 

Definition C.l A segment iy^,[a, b] of a stable manifold is the arc of the 

manifold between points a and h which includes these points. Correspondingly, 
Ws{aL, b) is the arc of the manifold between the points a and h excluding these 
points. 
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In the same manner one can define the segments iys[a, b) and iys(a, b]. The 
segments on the unstable manifold are defined similarly. 

Definition C.2 An initial segment of the branch of a stable manifold, W^{0, a], 
is the segment of which extends from the hyperbolic fixed point O to the 
point a G . 

Similarly the initial segment of the unstable manifold, VF^((9,a] can be de- 
fined. 

Definition C.3 A fundamental segment of a branch of a stable manifold, 
^^'''(a, T~^a], a G is the segment of between the point a and its 
T^^ -image, the point a is excluded and T~^a is included. Correspondingly, a 
fundamental segment of the unstable manifold, ^^^(a, Ta], is the segment of 
between the point a G (excluded) and its T-image (included). 

Some comments are in order. The maps T and T^^ transform fundamental 
segments into other fundamental segments. Namely, for any a G (or, 
correspondingly, a G W^) 

r{iy±(a,ra]} = W^±(Ta,T^a]; 
r{W^±(a,T-ia]} = iy±(Ta,a]. 

The whole branch of a stable (unstable) manifold is the union of all its fun- 
damental segments 

oo 

W^= U T"{iy„±(a,Ta]}, for any a G 

n=— oo 
oo 

py± _ y t''{W^{Tsl, a]}, for any a^eW^. 
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